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We consider a two-field model for inflation where the second order metric perturbations can 
be amplified by a parametric resonance during preheating. We demonstrate that there can arise 
a considerable enhancement of non-Gaussianity sourced by the local terms generated through the 
coupled perturbations. We argue that the non-Gaussianity parameter could be as large as /jvl ~ 50. 
Our results may provide a useful test of preheating in future CMB experiments. 



Introduction: Preheating, first realized in |lj and 
worked out in detail in j2, y, Bj] , may play an important 
role in understanding hot thermal Universe after the end 
of inflation. Preheating can occur during the coherent os- 
cillation of the homogeneous inflaton condensate if there 
exists a temporary vacuum instability due to which a 
resonant production of particles takes place. During this 
phase it is also possible to amplify the gravitational fluc- 
tuations or the metric fluctuation to the super Hubble 
scale 0, @ (for a general first order linear perturbation 
theory, see pj). So far the analysis has been limited to 
the first order metric and matter perturbations, although 
there has also been some attempts to understand the 
higher order effects /$). In some case the inflaton could 
also fragment into non-topological solitons, but then the 
inhomogeneities remain below the subhorizon scale [j| . 

Since CMB experiments such as WMAP have now 
reached a precision whereby cosmological models can be 
put to test, preheating may also move from the theoret- 
ical playground to observational scrutiny. In the present 
paper we propose that non-Gaussianities in the CMB 
fluctuations could provide a useful tool for studying pre- 
heating. 

The simplest single field inflationary models produce 
scale invariant Gaussian fluctuations, assuming that the 
initial state is the standard vacuum. However, non- 
Gaussianities are expected to be generated at some 
level, either because of inflaton self-couplings ^l), non- 
standard initial vacuum state |12( , or because of a host of 
other reasons often involvin g m odels beyond the simple 
single field inflation [iH. IbH Il5l] (for a review, see 
For example, in hybrid-type scenarios with two fields the 
fluctuations of the "waterfall" field may give rise to large 
non-Gaussianities during inflation [T^ . 

Major non-Gaussianities may also be generated after 
inflation Here we focus on non-Gaussianities dur- 

ing the preheating stage, although our approach is quite 
generic and applicable to any scalar condensate responsi- 
ble for reheating the Universe, such as a supersymmetric 
flat direction |18|. for which we may employ the formal- 
ism developed in |15|. We will consider a simple toy 
model of two fields, making some well-motivated physical 
assumptions regarding the motion of the coupled fields. 
Our starting point is an oscillating inflaton coupled to an- 



other field, which for simplicity is assumed to have a van- 
ishing vacuum expectation value. Therefore, classically 
the second field is sitting at the bottom of its potential. 
However, its fluctuation about the background solution is 
non- vanishing, giving rise to matter perturbations which 
naturally are higher order, mainly chi-squared fluctua- 
tions , which may source the metric perturbations and 
lead to non-Gaussianities. We illustrate that although 
the first order metric perturbations do not grow, for a 
long wavelength super-Hubble horizon modes the second 
order metric fluctuations can grow exponentially, high- 
lighting the potential for a large non-Gaussian signature 
arising from preheating. 

For simplicity, for obtaining the estimates we will ne- 
glect the expansion of the universe. This may seem rather 
restrictive, but since the homogeneous condensate oscil- 
lates coherently with a frequency larger than the Hubble 
expansion rate during the initial stages of preheating, the 
only relevant time scale is the mass of the coherently os- 
cillating field. Of course, for detailed numerical values, 
the expansion as well as backreaction effects will be im- 
portant. 

Basic Equations: Let us consider a two-field model 
with the potential 

t r 1 22,1222 

V = 2 m <^ + 2 9 V ° ' *• > 

where tp is the homogeneous scalar condensate coherently 
oscillating with a mass m^. We assume tp to be our 
inflaton. The a field will be created resonantly by virtue 
of the coupling term. 

For simplicity and for the sake of clarity we assume 
that the VEV of a vanishes, (a) = 0, which makes it pos- 
sible to obtain analytic approximations from the second 
order perturbation equations, as we show below. Such a 
situation occurs if the a field is driven to the minimum 
of its potential during inflation. 

The metric in our case is given by 

.900 = ~a(v) 2 (l + 2^ (1) + <P {2) ) , 
9oi = 0, 

9H = a(,7) 2 (l-2V^-^ 2) )%, (2) 
where we are using the generalized longitudinal gauge 
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and neglect the vector and tensor perturbations. Here r\ 
is the conformal time and a(rj) is the scale factor. We 
divide the fields into background and perturbations, 



<p = Mv) +5 {1) ip(r 1 ,x) + -S^<p(ri,x), 



(3) 



where the background value for a is assumed to vanish. 
The background equations of motion are found to be 



3?T = 



1 M 2 

1 / 2 i P 111 

2™2 , 



= <p^ + 2H<p' + a 2 m 2 (p 



(4) 
(5) 



while the <7o-equation is trivial, and Ji denotes the Hub- 
ble expansion rate in conformal time. We denoted the 
reduced Planck mass as M p ~ 2.4 x 10 18 GeV. 

The relevant 1st order perturbation equations can be 
written in the form Hal 



cW (1) + 2 (n 



<PoJ 

+2 (n'- ^j-n) ^ 



2.-2 S (i) a 



= 0. 



(6) 
(7) 



Note that there are no metric perturbations in Eq. JTJl. 
This is due to assuming a vanishing VEV for a. If it 
were non- vanishing, we would not be able to decouple 
the metric fluctuations from the perturbations in a field. 
Now the a part can be solved separately and for the rest 
the usual one field results apply. 



For our purposes it is sufficient to estimate the solution 
to be independent of k in the resonance band. There- 
fore, we estimate S (1) a,. lt = Ae~sq>(/j, B{t m lp t), where /i off = 
A*max/2 = q/i in the resonance band and £ (1> <T cff = oth- 
erwise; here A is an amplitude after the end of inflation. 

The amplitude depends on the effective mass of a 
field. If gipo < m tpi then the initial perturbations is 
given by A ~ S^a k ~ {H/V2~k^){k/aHf/ 2 - l/ where 
v = ^/9/4 — g 2 ip^/H 2 . For a narrow resonance regime, 
g < H/(p ~ H/M p <C 1, the spectrum follows that of 
the inflaton, where we assumed at the end of inflation 
ipo ~ M p . For a broad resonance regime, q > 1, and 
gtpo 3> m v , there the amplitude of the perturbations for 
a field will be suppressed compared to that of the infla- 
ton, see 0. 

Perturbations: The first order metric perturbation at 
the end of inflation for a single field is given by [j| , 



2M 2 



(10) 



where the quantities on the right-hand side are calculated 
at the horizon crossing. We obtained the final result by 
using slow roll equations of motion. Since g -C 1, then 
during inflation, we can effectively treat S^'cr ~ 
This determines the initial conditions before preheating 
starts. 

In the second order we are only interested in the grav- 
itational perturbation, whose equation can be written in 
an expanding background as 



2 ( W - ^-n 



8S l 



Jo 



local 



U a.i 



-local 



(ii) 



Parametric Resonance: Let us now turn to the para- 
metric resonance for 5 {1) a k . Throughout this section we 
neglect the expansion of the Universe, which results in 
a = 1 and n = t, i.e. conformal and cosmic time are 
equal, since the coherent oscillations begin only when 
m v 3> Ti. Hence the excitation of a k will follow the 
equation of motion valid for a narrow parametric reso- 
nance regime |3|, |5| , 



d 2 5 w 



0~k 



dz 2 



+ (A k - 2q cos2 z)S (1) a k = . 



(8) 



fe = k 2 jm 2 , +2q, q = g 2 <S> 2 /Am 2 v , and z = m v t. 



where A 

The amplitude of the oscillations is denoted by <F In a 
narrow resonance regime q < 1. In this regime the per- 
turbation grows exponentially as S (1) a k oc ex^>(p k m v t), 
where 



fi k = ^{q/2) 2 -{2k/m^-l) 2 



(9) 



Therefore, we have a resonance for the modes k when 
k- < k <k + with k± = (m lp /2)(l ± q/2). 



where the source terms J are quadratic combinations of 
first order perturbations; in particular, 



J^ C ^~^W) 2 + ^(S^) 2 



(12) 



Jo. non-local involves an inverse spatial Laplacian, thus 
rendering it non-local, while J7 rc st consists of metric and 
ip perturbations. Note that the left hand side of the above 
equation is identical to the first order equation Eq. ©. 

Fourier transforming J^iocai —> Jk (with our conven- 
tion f(x) = (l/27r) 3 J d 3 ke kx f(k)) we end up with the 
convolutions 



Jk 



M 2 (2n) 
1 5 2 V 1 



M 2 6a 2 (2tt) 3 



— Jd 3 k'6^a' h ,5^a' k _ k ' 

d 3 k'SVo- k ,SVo- k - k , . (13) 



The object is then to compute the convolutions. We are 
interested in their contributions at large scales and, to 
that end, we take the limit k — > 0, we also assumed 
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a = 1. Since the mode function 5 {1) a ott only depends on 
the magnitude of the vector k, the angular integration 
can be carried out trivially. The time derivative only 
produces a constant factor. Thus, we obtain 



Jk 



1 4tt 

a 2 6 2 V 4tt 



2fi 



2 2 
eff m f 



dk'k' 2 {5 {1) a' k ,) 2 

dk'k' 2 (S (1) a k ,) 2 
1 <5 2 Ul_ o 



M* 



M 2 5a 2 



in 



(14) 



dk' k' 



where in the last step we have assumed that the fc- 
dependence of the amplitude A can be ignored. If we 
are working in a narrow resonance regime with q < 1, 
the integral can be written as 



Jk- 



dk! k' 2 



1 fm v y 
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We can now write the source term as 



Jk 



47T 

(2^ 
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M 2 
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(15) 
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(16) 



where B = (q/8ir )(m ip /2) 3 A . It is worth noting that 
the source Jk , which we study at the large scale limit k ~ 
0, is actually generated by first order, local perturbations 
on much smaller scales (fc_ < k < fc + ). 

Non-Gaussianity: Consider Eq. (|lll) in fc-space. The 
homogeneous part is the same as in the first order. There- 
fore we know that the homogeneous solutions are well- 
behaved. Barring accidental cancellations, we may as- 
sume that the local terms we have considered are repre- 
sentative of the exponential behavior of the source, see 
Eq. I|12fl ■ There is also a non-exponential part which nat- 
urally becomes quickly insignificant. In order to estimate 
the behavior of <p {2) at large scales (k ~ 0) we neglect the 
expansion of the Universe and drop the terms with TL. 
The approximated metric perturbation then reads 



A: 



M 2 V 



l-± 
16 



+ [non-exponential source] , (17) 

Assuming that on average after some oscillations, the 
fraction <£o/Vo can be approximated by the frequency 
of the coherent oscillations ~ m Vl we readily obtain an 
exponential behavior for the solution of Eq. Q17[l. We 
may thus write 



2m 2 Lp {\- q/lQ) ■ B 
M 2 m 2 (l-q/4) ' 



qm v t/2 



(18) 



Let us use the following definition for the constant 

non-linearity parameter ffj L : (j> — <j> {1) + /tvlW 1 ') 2 - 
However, the definition is given in x-space and we 
have performed our calculations in fc-space. In princi- 

(2) 

pie we could transform qy, back to a;-space, but then 
we would need to know it for all k and we have eval- 
uated only the superhorizon mode k = 0. Instead, 
we can carry the definition of ffj L over to fc-space by 



,,(2) 



, where we have treated 



fx L as a constant. 



In the present scenario, where the 1st order pertur- 
bations are equivalent to that of a single-field case, we 
have S^a ~ S^ip right after inflation. Since, during 
inflation and preheating stays roughly constant, we 
immediately obtain an order of magnitude estimate from 
Eq. O 



4 2) 



,l- 9 /16 (M p 



NL 



<t>^)k 



l-?/4 V Ifo 



^ ) e N i' 2 ,(19) 



where we have written N = tu>, where N is the num- 
ber of oscillations during preheating and to is the fre- 
quency of the oscillations. On the average the frequency 
of the oscillations is given by to ~ m v . Inflation ends 
when ifo ~ M p , therefore the coefficient infront of the 
exponential is order one. The factor B in the coefficient 
of Eq. (|18fl . whose origin lies in the source terms, see 
Eqs. (|14I15I16|) cancels out completely from the contri- 
bution coming from * 4>^)k, because of the initial 
evolution for S^ip and are the same. The above 

expression should be compared with the observationally 

constrained one: Jnl — —/nl + (H/6)- 

The amplitude of the oscillations remains constant in 
our case. Therefore N is an indicative number, which 
is valid until the backreaction kicks in and shuts off the 
parametric resonance. Obviously N depends on a partic- 
ular potential and also on the expansion of the Universe. 
The parametric resonance would shift with the expansion 
as q ~ $ 2 t~ 2 and thus becomes narrower. However, for 
a simple single field chaotic type inflation model with 
V ~ (m 2 7 /2)(p 2 one would typically have many oscilla- 
tions within one Hubble time: uiH^ 1 ~ M p /$ 3> 1. It 
seems therefore that backreaction would be more deci- 
sive as far as the magnitude of the non-Gaussian ampli- 
tude fffij is concerned. This requires more study, but 
let us point out that in chaotic inflation backreaction 
becomes important after 10-30 oscillations Q. Hence 
N = O(10) might be a reasonable number, and should 
we for illustrative purposes choose q = 0.8, we would 
obtain f NL w e 4 w 55. This should certainly be at an 
observable level for the Planck Surveyor Mission. 

The expansion of the Universe changes the situation in 
two ways. First, the parametric resonance can be broad 
with q > 1 a time dependent quantity Second, be- 
cause of the the expansion the momenta and the oscil- 
lation amplitude redshift. However the amplitude of the 
first order metric perturbation still undergoes resonant 



4 



amplification as the momentum modes drift through the 
broad resonance regime p| . This ensures that the second 
order metric perturbations also grow exponentially, but 
one has to ensure that the amplitude of the initial per- 
turbations for a does not damp away during inflation. A 
detailed study would require numerical simulation, but 
nevertheless we may conclude that our result hints at 
the possibility of exciting the second order metric per- 



turbations during the first few oscillations of the infla- 
ton, hence linking preheating with possibly observable 
non-Gaussianities. 
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